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We investigate two types of chimera states, i.e., patterns consisting of coexisting spatially sepa¬ 
rated domains with coherent and incoherent dynamics, in ring networks of Stuart-Landau oscillators 
with symmetry-breaking coupling, under the influence of noise. Amplitude chimeras are charac¬ 
terized by temporally periodic dynamics throughout the whole network, but spatially incoherent 
behavior with respect to the amplitudes in a part of the system; they are long-living transients. 
Chimera death states generalize chimeras to stationary inhomogeneous patterns (oscillation death), 
which combine spatially coherent and incoherent domains. We analyze the impact of random per¬ 
turbations, addressing the question of robustness of chimera states in the presence of white noise. 
We further consider the effect of symmetries applied to random initial conditions. 

PACS numbers: 05.45.Xt,89.75.-k,05.40.Ca 


I. INTRODUCTION 

Chimera states are intriguing spatio-temporal patterns 
made up of spatially separated domains of synchronized 
(spatially coherent) and desynchronized (spatially inco¬ 
herent) behavior, although they arise in networks of com¬ 
pletely identical units. Kuramoto et al. discovered the 
chimeras in a network of phase oscillators with a sim¬ 
ple symmetric non-local coupling scheme [T]. Within 
these “classical” chimera states, coherent domains of pe¬ 
riodic in-phase oscillations coexist with incoherent do¬ 
mains, characterized by a chaotic behavior in time. In 
Greek mythology, the chimera is a fabulous creature with 
a lion’s, a goat’s, and a snake’s head. As the counter intu¬ 
itive dynamical state is as well composed of incongruous 
parts, it was named after the beast HIE]- Kuramoto’s 
finding initiated a broad wave of theoretical and empiri¬ 
cal investigations HE?]. The first experimental evidence 
on chimera states was presented only one decade after 
their theoretical discovery [28ti38] . In real-world sys¬ 
tems chimera states might play a role, e.g., in the uni- 
hemispheric sleep of birds and dolphins [39] , in epileptic 
seizures m, in power grids [41] , or in social systems [42] . 

The basin of attraction of the chimera states is 
typically relatively small compared to that of the 
synchronized state. This explains why they were not 
detected for a long time. Further, it is the reason why for 
many investigations specially prepared initial conditions 
are used. An interesting feature of the chimera states 
is that they are often long living transients towards the 
in-phase synchronized oscillatory state. By this, the 
coupling between the oscillators introduces a time scale 
much larger than the oscillation periods of each single 
oscillator. It was theoretically predicted [13] and experi¬ 
mentally confirmed [37] that the lifetime of the chimera 


states of phase oscillators grows exponentially with the 
system size. This illustrates impressively that they are 
not simply a temporary trace of the initial conditions, 
but are a persisting phenomenon. Recently, the inves¬ 
tigation of chimera states has been generalized towards 
networks of elements which have more complicated local 
dynamics 0 EH EE] • In particular, nodes which not 
only involve phase but also amplitude dynamics are 
considered. As described in [20] [21], in such systems 
amplitude-mediated chimeras can be found, which show 
a chimera behavior with respect to the phases as well 
as with respect to the amplitudes. In this work, we 
will consider a different type of chimera states. They 
are characterized by strictly correlated phase dynamics 
throughout the whole network, but coexisting domains 
of coherent and incoherent amplitude dynamics. These 
amplitude ehimeras were first described in [43] . A crucial 
difference to classical phase chimeras is that the spatial 
incoherence does not imply chaotic behavior in time. In 
fact, all nodes of an amplitude chimera perform periodic 
oscillations, but in the incoherent domain the spatial 
sequence of the positions of the centers of oscillation 
is completely random. It has been noted [44] that 
amplitude chimeras are long living transients, however, 
the dependence of their lifetimes upon the system size is 
not known. This will be one of the questions addressed 
in this work. 

It has been shown that the occurrence of the amplitude 
chimeras in a network of Stuart-Landau oscillator is 
restricted to the case of non-local coupling [43] m]. 
Furthermore, the breaking of the rotational symme¬ 
try of the single oscillator is crucial. In this work, this 
symmetry breaking is realized through a coupling term 
which only involves the real parts instead of the complex 
variable 2 :, in accordance with [43] . 

This symmetry-breaking coupling introduces a nontrivial 
spatially inhomogeneous fixed point (steady state) into 
the system [46], which is typically associated with a 
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pitchfork bifurcation. As a result, oscillation death 
patterns occur in addition to the oscillatory states. 
Oscillation death is a phenomenon appearing in oscil¬ 
lator networks, when the coupling induces a splitting 
of a (trivial) homogeneous fixed point into at least two 
distinct branches gillTl. This contrasts with amplitude 
death, where the quenching of the oscillations is caused 
by stabilization of an unstable spatially homogeneous 
fixed point which is already present in the uncoupled 
case [48]. When the newly born inhomogeneous fixed 
point is stabilized, typically through an inverse Hopf 
bifurcation, the oscillatory dynamics can be suppressed. 
In general, however, oscillation death states can coexist 
with oscillatory solutions. Amplitude chimeras and 
oscillation death are two different phenomena related 
to symmetry breaking that occur in the Stuart-Landau 
oscillator network. Interestingly, a novel type of steady 
state coherence-incoherence patterns combining features 
of both phenomena, called chimera deaths has been 
discovered in this system [43]. Like amplitude chimeras, 
these patterns are characterized by coexisting domains 
of coherent and incoherent behavior in space. They have 
recently been extended to globally coupled networks 
[49] ED]. 

An important question is whether transient chimera 
patterns persist under perturbations. The robustness of 
classical phase chimeras with respect to heterogeneous 
natural frequencies of the oscillators [SIEIIISS or hetero¬ 
geneous connections [53] was investigated. It was found 
that the specific way in which links are removed from 
the network defines how strongly the chimera lifetimes 
are affected. Also in the FitzHugh-Nagumo model, the 
robustness of chimeras with respect to heterogeneous fre¬ 
quencies and topologies was studied [25] . Apart from the 
present work, the robustness of chimera states with re¬ 
spect to external noise has so far not been considered. 
This problem is especially relevant in the light of experi¬ 
mental realizations where noise is inevitable. It is known 
that even at a relatively low intensity, noise can signif¬ 
icantly influence the behavior of a nonlinear dynamical 
system. Under certain conditions noise can cause the in¬ 
crease of coherence, e.g., in coherence resonance, and on 
the contrary it can also induce irregular, chaotic dynam¬ 
ics. Stochasticity appears due to the intrinsic fluctua¬ 
tions in a system, or alternatively can be implemented as 
an external random control force. The influence of noise 
strongly depends on its characteristics such as intensity, 
for example, which can be treated as bifurcation param¬ 
eters of a stochastic system. Therefore, the investigation 
of noisy dynamics is on the one hand significant for the 
understanding of the processes occurring in nature and 
on the other hand relevant from the point of view of con¬ 
trol. 

In this work, we investigate how external noise 
influences the occurrence of chimera death states and 
amplitude chimeras, and how it affects their properties 
and lifetimes. For this purpose, we introduce novel order 


parameters to characterize chimera states, and discuss 
the impact of initial conditions. In order to explain our 
numerical findings, we propose that amplitude chimeras 
represent a saddle-state in the underlying phase space. 

The paper is organized as follows. In Sect. II we in¬ 
troduce the model. In Sect. Ill we analyze amplitude 
chimeras and chimera death without noise and introduce 
a global measure which can distinguish between chimera 
patterns and completely coherent dynamics like in-phase 
synchronization or traveling waves and can be used to 
quantify the lifetime of chimeras. Further, we focus on 
the role of the initial conditions and the dependence of 
the chimeras upon coupling parameters and system size. 
In Sect. IV the influence of noise upon the transient 
times of chimeras is investigated. Sect. V summarizes 
the results in terms of the different regimes of the cou¬ 
pling parameters, and Sect. VI draws conclusions. 


II. MODEL 

We consider a network of N Stuart-Landau oscillators 
[H |43l im (46] [54] under the impact of external white 
noise ^j{t). The local deterministic dynamics of each 
node j G {!,..., V} is given by Zj = f{zj)^ with the 
normal form of a supercritical Hopf bifurcation 

f{zj) = {X + iuj-\zj\^)zj, (1) 

where Zj = Xj -\- i i/j =rje^^^ G C, with Xj^yj^Vj, <pj G M, 
and A,cc; > 0. At A = 0 a Hopf bifurcation occurs, so that 
for A > 0 the single Stuart-Landau oscillator exhibits 
self-sustained oscillations with frequency uj and radius 
Tj = a/A, and the unique fixed point Xj =0, yj = 0 is 
unstable. 

We investigate a ring of N non-locally coupled Stuart- 
Landau oscillators, where each node is coupled to its P 
nearest neighbors in both directions with the strength 
cr > 0, and is subject to noise of intensity D > 0 : 

j+P 

= fizj) + ^ V (2) 

k=j-P 

where j = 1, 2,..., V and all indices are modulo N. The 
coupling and the noise are only applied to the real parts. 
In many real-world oscillator systems only one variable 
is accessible for coupling. For instance, the diffusive cou¬ 
pling through only one variable naturally arises in bio¬ 
logical and biochemical systems. The simplest coupling 
scheme is hence to choose the real parts of the com¬ 
plex variable z. Here ^j{t) G M is additive Gaussian 
white noise [55], i.e., {^j{t)) =0, Vj, and {^i{t)^j{t')) = 
6ij6{t — t'), Vi, j, where 6ij denotes the Kronecker-Delta 
and S{t — t') denotes the Delta-distribution. Hence the 
noise is spatially uncorrelated. 
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FIG. 1. (Color online) Amplitude chimera state in system 0 
with N = 200 nodes, for coupling range P/N = 0.04 and cou¬ 
pling strength cr = 18: (a) snapshot (top: Xj, bottom: yj), (b) 
space-time plot, (c) phase plot in the complex plane: trajec¬ 
tories of 12 nodes of the coherent domains (unit circle, blue) 
and 12 nodes of the incoherent domains (red and green), the 
arrows indicate the direction of the motion. Initial condition: 
see section Other parameters: F> = 0, A=l, Ct; = 2. 


III. DETERMINISTIC CHIMERA PATTERNS 
AND THE ROLE OF INITIAL CONDITIONS 


As shown in [43] in the absence of noise, various differ¬ 
ent states can be found in the network given by Eq. 
Which particular state actually arises, depends on the 
specific values of the coupling parameters and the ini¬ 
tial conditions, as Eq. <§ describes a multistable system. 
Among the possible states, two different types of asymp¬ 
totically stable states can be found, on the one hand oscil¬ 
latory states, and on the other hand steady state patterns 
which are related to oscillation death. The latter are rep¬ 
resented by completely coherent or completely incoher¬ 
ent oscillation death patterns, as well as by chimera death 
patterns. The asymptotically stable oscillatory states ap¬ 
pear in two different spatio-temporal patterns: in-phase 
synchronized oscillations and coherent traveling waves. 
Besides these, long lasting oscillatory transients with in¬ 
teresting features occur, i.e., amplitude chimera states. 
In this work we demonstrate that all these states can 
also be observed under the influence of noise. 


A. Amplitude chimeras and chimera death 
configurat ions 

Before an asymptotic oscillatory state (a completely 
in-phase synchronized oscillation or a coherent travelling 
wave) is approached, amplitude chimera states can ap¬ 
pear as long transients, potentially lasting for hundreds 
or even thousands of oscillation periods. In contrast to 
classical chimeras, during the chimera all nodes (includ¬ 
ing the ones within the incoherent domains) oscillate with 
the same approximate period, T = ^, and a spatially 



FIG. 2. (Golor online) Snapshots of coherent oscillation death 
states: (a) coupling strength cr = 18 and coupling range 
P/N = 0.14 (2-cluster), (b) a = 8, P/iV = 0.04 (10-cluster). 
Initial condition: nodes 0 < j < 24 and 50 < j < 74 
are set to {xj,yj) = (0.1,—1), all other nodes are set to 
{xj,yj) = (—0.1,+!). Other parameters: N = 100, P = 0, 
A = l, a; = 2. 


correlated phase, but they show spatially incoherent be¬ 
havior with respect to the amplitudes in part of the sys¬ 
tem. Figure shows an exemplary amplitude chimera 
configuration. The nodes within the two coherent do¬ 
mains (here 13 < j < 85 and 113 < j < 185) perform 
synchronized oscillations, all with the same amplitudes. 
The coherent domains always appear pairwise, such that 
for every time t, all nodes within one coherent domain 
have opposite phases than all nodes of the other, antipo¬ 
dal domain. Hence they always fulfill the “anti-phase 
partner” condition Zj = —5 j A, assuming 

even N. As visible in Eig.lH (c), the trajectories in the 
complex plane of all nodes are cycles, illustrating that 
all nodes have periodic dynamics in time. This is a fun¬ 
damental difference between the classical chimera states 
where a part of the network demonstrates chaotic tem¬ 
poral behavior. The nodes of the coherent domains all 
oscillate on a perfect circle around the origin. Both co¬ 
herent domains are represented by one single blue line 
in Eig.[^ (c), which is at the same time the trajectory 
of all nodes when the completely in-phase synchronized 
oscillatory solution is approached. The two antipodal 
coherent domains are separated by incoherent domains. 
There, neighboring nodes can be in completely different 
states at a given time t. Their trajectories are deformed 
circles, whose centers are shifted from the origin. In the 
incoherent domain the sequence of nodes that oscillate 
around centers in the upper and lower half-plane, respec¬ 
tively, is completely random. The ratio of the domain 
sizes can vary strongly depending upon the parameters 
as discussed below, but does not depend upon the initial 
conditions. Transient amplitude chimeras with very nar¬ 
row incoherent domains can be observed, as well as with 
broad ones. 

If the coupling strength and coupling range exceeds cer¬ 
tain values, the oscillations of the Stuart-Landau nodes 
can be suppressed due to the stabilization of a new in¬ 
homogeneous steady state created by the coupling. In¬ 
stead of performing oscillations, each node approaches a 
fixed point close to one of the following two branches: 
(x*^,^*^) ^ (—0.1,4-0.85) or (x*^,^*^) « (-j-O.l, —0.85) 
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FIG. 3. (Color online) (a),(b) 1-cluster chimera death (1-CD) 
for coupling range P/N = 0A, (c),(d) 3-cluster chimera death 
(3-CD) for P/N = 0.2. Snapshots Xj,yj are shown in panels 
(a),(c), space-time plots in panels (b),(d). Initial condition: 
see section [ill C| Parameters: A/" = 200, a = 18, i4 = 0, A=l, 

ct; = 2. 


for all other configurations. However, this order parame¬ 
ter cannot be usefully applied to determine the transient 
time of an amplitude chimera. We note that the in-phase 
synchronized oscillation is not the only asymptotically 
stable solution the system can approach; traveling waves 
are also possible, but R{t) cannot distinguish between 
amplitude chimeras and coherent traveling waves since 
both give R{t) « 0. 

In order to define an appropriate order parameter, we in¬ 
troduce the center of mass of each oscillator 

jm- 


1 


27r 


vr- = fj 

f yj{t) dt : with T - 

5 

UJ 

(4) 

and an analogous 

definition of . The 

shift 

of the 

center of mass from the origin is given by: 



r^c.m. 

3 

II 

to 

+ 

to 


(5) 


(for A = 1), and remains there for all times. The oscil¬ 
lation death states exhibit a huge variety of spatial pat¬ 
terns, including multiple coherent and multiple incoher¬ 
ent oscillation death states [43l[44]. Two exemplary con¬ 
figurations of completely coherent oscillation death pat¬ 
terns are shown in Fig.[^ (2-clnster and 10-cluster oscil¬ 
lation death), see also [47]. The oscillation death regime 
is characterized by very high multistability. Among the 
oscillation death states, chimera death patterns can be 
found, which combine the characteristics of both phe¬ 
nomena: chimera state and oscillation death. These pat¬ 
terns consist of coexisting domains of coherent and in¬ 
coherent populations of the inhomogeneous steady state 
branches. Within the incoherent domains, the popula¬ 
tion of the two branches (upper and lower) follows a ran¬ 
dom sequence, as for example visible in Fig.|^ Within the 
coherent domains, the number of clusters of neighboring 
nodes that populate the same branch of the inhomoge¬ 
neous steady state can vary. An m-cluster chimera death 
state (m-CD), with m G {1, 3, 5, 7, 9,...}, is characterized 
by the occurrence of m clusters within each coherent do¬ 
main. The coherent domains always appear pairwise with 
anti-phase symmetry Zj = similarly to the co¬ 

herent domains of the amplitude chimera configurations. 
Our numerical results confirm that the stable oscillation 
death patterns fulfill the “anti-phase partner” condition. 


B. Order parameter 


The global Kuramoto order parameter mm 


m 


1 

N 


N 

i=l 


(3) 


can be used to detect if the system is in the in-phase syn¬ 
chronized oscillatory state, where R{t) = 1. It is smaller 


These measures allow to distinguish between coherent 
and incoherent domains: and hence also 

vanish for nodes within the coherent domains of 
the amplitude chimera. In contrast, all nodes within 
incoherent domains have finite values. Here, the spatial 
sequence of positive and negative signs of y^-'^- is 
completely random. Therefore, these quantities can 
serve as local order parameters. 

We propose the global order parameter based on their 
second moments 


C2{t) = 


\ 


1 ^ 


(6) 


i=l 


to detect the transition from an amplitude chimera to a 
coherent oscillating state. Figure [^ shows two examples 
with different initial conditions for such a transition, one 
where the approached stable asymptotic solution is an 
in-phase synchronized oscillation, and one where it is a 
coherent traveling wave. Both figures depict a space-time 
plot (upper panel), and the temporal evolution of the 
global order parameters R{t) and C 2 (t) for a time series 
around the transition. As is clearly visible, in contrast to 
the global Kuramoto order parameter, C 2 is capable to de¬ 
tect both types of transition. Another advantage of this 
order parameter is that its value solely depends on the 
relative size of the incoherent domains of the amplitude 
chimera, which in turn depends only on the system pa¬ 
rameters. In particular, it does not depend on the pattern 
of disorder within the incoherent domains, and hence not 
on the specific realization of the initial condition. This 
means that the value C 2 is a characteristic quantity of a 
chimera, which distinguishes it from other patterns. We 
define the relative number of nodes within the incoher¬ 
ent domains as the total number of nodes within both 
incoherent domains normalized by the overall number of 
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FIG. 4. (Color online) Transition from an amplitude chimera state to (a) in-phase synchronized oscillation, (b) coherent 
traveling wave. Upper panels: space-time plots, lower panels: time series of the order parameters: global Kuramoto order 
parameter R{t), second moments C 2 , C 2 with a = 0.0361. Parameters iV = 100, P/N = 0.04, a = 14, D = 0, A = 1, a; = 2. In (a) 
and (b) different realizations of the initial conditions shown in Fig. 5(b) with Q = N/4 are used. 




FIG. 5. (Color online) Specially prepared initial conditions 
for amplitude chimeras: (a) point-symmetric type, (b) fully- 
symmetric type. Top panels: Xj, bottom panels: yj. Dashed 
lines indicate the point-symmetry about the center, vertical 
solid lines and arrows (green) indicate the axial symmetry 
within both network halves. System size: N = 100. 


nodes within the whole network. To calculate the rela¬ 
tive number of nodes within the incoherent domains, we 
propose a modified version: 

1 ^ 

= ( 7 ) 

i=i 

with the Heaviside function 0(x), and some fixed 
threshold a > 0. The Heaviside function is non-zero 
only for those nodes where is above this 

threshold. For the choice a = 0.0361, the predicted 
number agrees very well with visual determination of 
the incoherent domain size, which we have checked with 
numerous amplitude chimera configurations in a wide 
range of system parameters. 


C. Role of initial conditions 

Eq.§ is known to describe a multistable system [43] . 
Both types of chimera states appear in coupling pa¬ 
rameter regimes, where other oscillation death patterns 


and coherent oscillatory states can be found as well. 
In order to increase the probability of finding chimera 
states, we use specially prepared initial conditions. A 
very simple initial condition that produces transient 
amplitude chimeras in a certain parameter regime (of 
about 0.01 < P/N < 0.05, a < 33), is when all nodes of 
one half of the network (1 < j < y) set to the same 
value {xj,yj) = (xo,^o) (excluding the choice (0,0)), 
and the rest is set to (xo,^o) ~ Hence, 

amplitude chimera states can evolve out of an initial con¬ 
figuration that only consists of two completely coherent 
parts. We choose the values (xj, y^) = (\/0^, — \/^), so 
that the nodes start on the limit cycle with r = \f\ = 1, 
which is the solution for the in-phase synchronized 
oscillation. The amplitude chimera lifetime nevertheless 
appears to be of the same order for other values (e. g. 


By adding random numbers to yj, we construct a 
more general class of specially prepared random initial 
conditions for amplitude chimeras. In particular, we add 
a random number drawn from a Gaussian distribution 
with variance V to yj of the Q nodes on the left and 
on the right side of the borders between both halves (at 
j = Y = indicated in Fig.|^ (a), with Q e N 

and 0 < Q < Besides the range Q of incoherence, we 
also vary U > 0. For a proper choice of the two initial 
condition parameters (Q and V), we obtain amplitude 
chimeras. Using the achieved amplitude chimera lifetime 
as a quality measure for the initial condition, we compare 
multiple realizations of the specially prepared random 
initial conditions for the deterministic system with 
P/N = 0.04 and N = 100. We observe that among all 
considered kinds of initial conditions (different choices 
for Q and U, symmetry conditions, xj randomized as 
well, a different underlying distribution for the random 
numbers), the applied symmetry of the initial condition 
has the greatest effect upon the transient time. 
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FIG. 6. (Color online) Transient times of amplitude chimera 
states ttr for 20 realizations of specially prepared random ini¬ 
tial conditions (no symmetry, point-symmetry, full symmetry, 
depicted by different symbols) with 10, V — 2. Horizontal 
solid lines: mean values. Dashed line: Ur for the initial con¬ 
dition with y = 0 (no randomization). System parameters: 
iV=100, cr = 14, P/iV = 0.04, D = 0, A = l, u = 2. 


Figure shows the transient times and their mean 
value (solid lines) for multiple realizations of the initial 
conditions following three different symmetry schemes. 
For the particular choice (Q = 10, V = 2) all symmetry 
types lead on average to shorter lifetimes than an initial 
condition without random component (black dashed 
line). For the initial configurations without symmetry, a 
random number is chosen independently for each node 
within the four incoherent intervals. These configu¬ 
rations clearly create the shortest amplitude chimera 
lifetimes, lasting at most for a couple of oscillation 
periods. This symmetry type also leads to the shortest 
transients in other regimes of Q, and V (not shown 
here). In contrast, for the point-symmetric initial 
conditions, we mirror the random numbers used for 
j e {1,..., y} with respect to the center j = 0, yj = 0, 
and use their negative counterparts for the second half. 
We hence only generate 2Q random numbers in total. 
The initial configurations are point-symmetric with 
respect to the center, see Fig.[^ (a). The lifetimes of 
the occurring amplitude chimeras are much longer than 
in the non-symmetric case. However, the symmetry 
type which leads to the longest lifetimes is the one 
referred to as full symmetry; the initial conditions fulfill 
two symmetries: The randomly chosen values of the 
positions of the nodes within the first incoherent interval 
1 < J < Q are mirrored to the nodes ^ — Q < j < y, 
by setting: Zj = znj^^_-. To obtain the positions of the 
second network half, then a phase shift of tt is applied, 
such that the “anti-phase partner” condition is fulfilled 
{zj = and j mod N). Thus, we only generate Q 

different random numbers in total. The configurations 
are again point-symmetric with respect to the center, 
and have an additional axial symmetry with orthogonal 
axes through j = ^ and j = ^, as indicated in Fig.|^ 
(b). Of course, the simple initial condition with no 
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FIG. 7. (Golor online) (a) Fully-symmetric random initial 
condition with range of incoherence Q = "^ and variance V = 
0.1. (b) Transient times of amplitude chimeras Ur for 40 

realizations of the initial condition shown in (a): circles: V = 
0.5, triangles: V = 0.1. Horizontal solid lines: mean values, 
dashed line: Ur for H = 0. System parameters: N = 100, 
cr = 14, P/iV = 0.04, D = 0, A = l, uj = 2. 


randomization also fulfills these symmetry conditions 
and can therefore be regarded as one special type of 
the fully-symmetric specially prepared initial conditions 
(with y = 0). We have further tested another type of 
initial condition that solely fulfills the anti-phase partner 
condition: zj = —Zj^n^ but has no other symmetries. 
This type of initial condition also certainly leads to 
transient amplitude chimeras, but only within very 
narrow ranges of Q and V. For (5 = 10, V = 2^ the mean 
lifetime (of about Ur ~ 49), is only slightly increased 
compared to the non-symmetric initial condition (not 
shown here). 

Since the symmetry which is applied to the initial 
conditions remains preserved during the dynamic evo¬ 
lution, this observation means that the fully-symmetric 
amplitude chimeras are most stable and have the longest 
lifetimes. 

By decreasing the variance in the interval 0.1 < ^ < 2, 
the mean amplitude chimera lifetimes increase. In the 
range of small variances of about V < 0.5, amplitude 
chimeras occur for all choices of the incoherence range 
(5, and the particular choice of Q does not infiuence the 
transient time much. For (5 = ^, all nodes are random¬ 
ized, see Fig.[^(a), which appears to be a natural choice. 
Fig.j^b) shows the corresponding transient times belong¬ 
ing to a set of 40 realizations of the specially prepared 
random initial condition with Q = ^ and V = 0.5, and 
for a set with V = 0.1. The mean transient times are 
much longer than for 1/ = 2 (cf. Fig.[^. They are at least 
of the same order (and can be larger) as the transient 
time for the simple initial condition with no randomiza¬ 
tion, V = 0 (dashed black line). For the choice V = 0.1, 
the transient times are increased as compared to 17 = 0.5. 
We use the fully-symmetric initial conditions with Q=^ 
and 17 = 0.1 for all investigations presented in this paper. 

Besides oscillatory states, oscillation death states can oc¬ 
cur in a large variety of different spatial patterns. Our 
numerical results suggest that in the appropriate parame- 
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P 

FIG. 8. (Color online) Relative size of the incoherent domains 
in the plane of coupling strength <7 and nearest neighbors 
number P. Initial condition: a fixed snapshot of an amplitude 
chimera. The light grey area denotes chimera death (CD). 
Other parameters: A^=100, D = 0, A=l, = 


ter regime every amplitude chimera snapshot can be used 
as initial condition to certainly produce a chimera death 
state. How many clusters in the coherent domain of the 
chimera death pattern occur, depends on the initial con¬ 
dition as well as on the parameter choice (see section [V|). 


D. Relative size of the incoherent domains 

The ratio of coherent and incoherent domains of an 
amplitude chimera can be measured by the number of 
nodes within the incoherent domains divided by the total 
number of nodes. It can be calculated using the order 
parameter C 2 (t). Interestingly, in the time evolution of 
the system, this quantity typically fluctuates by less than 
one order of magnitude, during the chimera, as well as 
during the coherent synchronized oscillation, or during 
the traveling wave. This becomes visible by considering 
the temporal evolution of C 2 (t), as exemplarily shown in 
Fig. El In contrast to R(t) and C 2 (t), the values of C 2 (t) 
show a sharp drop to zero when the amplitude chimera 
disappears, illustrating that the relative size of the 
incoherent domains remains constant until the last oscil¬ 
lation period before the transition to complete coherence. 

Note that the relative size of the incoherent domains 
is not related to the choice of Q, but is completely inde¬ 
pendent of the initial condition. Even if the incoherent 
domain is chosen larger or smaller initially, it evolves 
into a domain with a characteristic size corresponding to 
the given set of parameters. This is an important general 
distinguishing factor of chimeras which emphasizes that 
they are self-organized patterns. 


(f) 



FIG. 9. (Color online) Relative size of the incoherent domains 
of amplitude chimeras vs. coupling strength cr, measured by 
C2(t=100). Squares, triangles, and disks correspond to three 
different initial conditions. Solid line: linear fit to disks. Sys¬ 
tem parameters: iV=100, P/iV = 0.04, T) = 0, A=l, u; = 2. 


In contrast, the relative size of the incoherent domains 
does depend on the choice of the coupling parameters. 
This is shown in Fig.for coupling parameters varied in 
the range 1 < P < 20, and 5 < a < 30. For each choice or 
the coupling parameters, the same snapshot of an ampli¬ 
tude chimera is used as initial condition. The maximum 
simulation time is t = 5000. The relative size of the in¬ 
coherent domains varies over a wide range from about 
0.14 to about 0.6, i.e., we observe amplitude chimeras 
where just 14% of the nodes show a spatially incoher¬ 
ent behavior, up to configurations with 60%. At the 
border between the oscillatory regime and the chimera 
death regime, we detect the amplitude chimeras with the 
largest incoherent domains. The width of the incoher¬ 
ent domains appears to increase linearly with P, as well 
as with cr. In Fig.[^ the relative size of the incoherent 
domains in dependence upon the coupling strength is de¬ 
picted for three different initial conditions. The actual 
spread of the values is due to the fact that, firstly, C 2 (t) 
can assume only multiples of and, secondly, the rela¬ 
tive number of nodes can change only in multiples of 
due to the symmetry in the initial conditions (see section 


IIIC). As opposed to the linear increase with cr and P, 


the relative size of the incoherent domains is roughly con¬ 
stant for all system sizes. For an exemplary choice of the 
coupling parameters. Fig. 10 shows (a) snapshots and (b) 
the corresponding shifts from the origin of the centers of 
mass, for five different system sizes N. 


E. System size dependence 

For classical phase chimeras it has been established 
that the transient nature of the chimera states is a 
finite-size effect [ails]. In the limit of an infinitely large 
system, the phase chimera is stable. This is opposite to 
the system size dependence of the amplitude chimera 
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FIG. 10. (Color online) Snapshots of amplitude chimeras for 
different system sizes iV, but constant coupling range P/N = 
0.04: (a) yj (b) shift of the center of mass from the 

origin. Parameters: t = 140, cr = 14, D = 0, A = 1, uj — 2. 
Initial conditions: see section [Til Cl 


FIG. 11. (Color online) Transient times of amplitude chimeras 
itr VS. system size A/", for the constant coupling range PjN — 
0.04. Disks: medians of itr over 20 initial conditions (IC), 
solid line: fit to Eq. Q with ai 3802.29, and a 2 ~ 596.23. 
Triangles: Ur for a single initial condition. Other parameters: 
cr=14, D = 0, A=l, uj = 2. 


by a sweeping front which travels through the whole sys¬ 
tem. Rather, the incoherent oscillators of the amplitude 
chimeras catch up locally with the coherent ones to syn¬ 
chronize, which is independent of the global system size 
as long as this size is large enough [44]. This reasoning 
is in line with our observations in Sect. mm that the 
relative size of the incoherent domains firstly, is roughly 
constant for all N (Fig.[TQ|), and secondly, stays constant 
during the temporal evolution of a chimera until the last 
period before the transition to the asymptotic solution 
and then drops immediately to zero (Fig.|^. 


states, where we find that their lifetime decreases with N. 


IV. INFLUENCE OF NOISE ON TRANSIENT 
TIMES 


For various system sizes V, Fig. [TT] shows the median 
of the transient times for a constant coupling strength a = 
14 and a constant coupling range P/V = 0.04, averaged 
over 20 initial conditions (black disks). The data for an 
additional single initial condition are also shown by cyan 
triangles. Both data sets can be well described by the 
following dependence of the amplitude chimera lifetime 
on the system size: 


t,AN) = 


0^1 V , 

J + 0 ^ 2 - 


(8) 


The black curve in Fig. is a fit with ai 3802.29 
and 0^2^596.23. This indicates that amplitude chimeras 
should exist, albeit with shorter lifetime, in systems of ev¬ 
ery size. That means that they would also appear in the 
“thermodynamic limit” of infinitely large systems with 
a finite lifetime (for instance, Ur ~ 600 for V ^ oo, for 
the type of initial condition used in Fig. 11). These re¬ 


sults show that the mechanism of the transient is differ¬ 
ent from phase chimeras, where the transition is triggered 


In this section we will study the robustness of chimera 
states with respect to external noise. By using the same 
initial conditions which lead to amplitude chimera states 
and chimera death in the deterministic case, we also ob¬ 
serve these states in Eq. ([^ in the presence of noise in 
a wide range of the coupling parameters. The stochas¬ 
tic equations were integrated with the well-established 
Euler-Maruyama scheme with stepsize 10“^, see p. 340 
in [45|. Figure [T^ shows one exemplary configuration for 
an amplitude chimera which occurs in a system under 
the impact of noise of intensity P = 5 • 10“^. 


In general, the transient times of amplitude chimeras 
decrease with increasing noise intensity. Figure [T^ shows 
the average transient times and the corresponding stan¬ 
dard deviations in dependence of the noise intensity P, 
for three choices of the coupling strength cr, in a semi- 
logarithmic plot. The average is over 50 different fully 
symmetric initial conditions (with Q = ^, see section 


III C) drawn from different realizations of the associated 
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node j 

FIG. 12. (Color online) Amplitude chimera for noise intensity 
F> = 5 • 10“^: (a) snapshot (top: Xj^ bottom: (b) space- 

time plot. Parameters: N = 200, P/N — 0.04, a = 19, A = 1, 
uj = 2. Initial condition: see section mn 



FIG. 13. (Color online) Transient times of amplitude chimeras 
ttr VS. noise strength D (log-scaled) for different values of cou¬ 
pling strength a. Symbols: average over 50 fully symmetric 
initial conditions (with Q = A ^ each associated with a differ¬ 
ent realization of the random force ^(t); error bars: standard 
deviations; lines: linear fits from Eq. (1^. Inset: slope vs a. 
Parameters: iV=100, P/N = 9.94:, A=l, uj = 2. 



FIG. 14. (Color online) Transient times of amplitude chimeras 
Ur averaged over 50 initial conditions and noise realizations 
(the same set of initial conditions as used in Fig. 13): (a) Ur 
vs. coupling strength a for different noise intensities. Sym¬ 
bols: mean transient times; error bars: standard deviations. 
The lines serve as a guide to the eye. (b) Ur in the plane of 
coupling strength a and noise intensity D. Other parameters: 
iV=100, P/N = 9.94, X=l, uj = 2. 



FIG. 15. (Color online) Schematic phase-space structure of an 
amplitude chimera (AC) as a saddle-point. Thick solid lines: 
stable directions, thick dashed lines: unstable directions, thin 
solid lines: different trajectories, with arrows denoting the 
direction of time evolution. Grey shaded region: scheme of 
amplitude chimera configuration, green (dark grey) disk: set 
of initial conditions (IC), yellow (light grey) area: impact of 
Gaussian white noise. 


random distribution. For each one of those realizations of 
the initial conditions, a different realization of the Gaus¬ 
sian white noise ^j{t) is considered. The average tran¬ 
sient times show a clear linear decrease as a function of 
the logarithmic noise intensity. This behavior is found 
throughout the range 6 < cr < 24, i.e., Ur = —^ ln{D)4-r] 
with slope — ^ and axis intercept r]. This gives the scaling 
law 

(9) 

The lines in Fig.[^ show the linear fits, and the inset 
depicts the slope in dependence on the coupling strength 
a. For the same set of 50 initial conditions, Fig.p^a) 
depicts the mean transient time in dependence of the 
coupling strength for four different noise intensities 
D, and Fig.p^b) shows a color-coded density plot of 


the mean transient times of amplitude chimeras in the 
{a, D)-plane. The transient times generally decrease 
with increasing noise, and increase with increasing 
coupling strength up to a saturation value at about 
(T'^lb. The vertical error bars in panel (a) show that the 
transient times are less sensitive to the initial condition, 
the larger the noise is. We generally find that the spread 
of the amplitude chimera lifetimes for different initial 
conditions (and different noise realizations), is smaller 
with increasing noise strength. 

Transient amplitude chimeras can last for thousands 
of oscillation periods until they disappear. Even under 
disturbance by external noise they persist for a sig¬ 
nificant time. Noise does not essentially change their 
spatial configuration. If noise throws the system onto an 
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adjacent trajectory in the underlying high-dimensional 
phase space of the network, this does not normally lead 
to a flow into a completely different direction in phase 
space. Geometrically speaking, this shows that there are 
some attracting directions in phase space along which 
the system dynamics is pushed towards the amplitude 
chimera. Furthermore, amplitude chimeras can evolve 
out of initial configurations that do not show the charac¬ 
teristic coexistence of coherent and incoherent domains 
(see paragraph III). In fact, they can be observed when 
completely incoherent initial configurations are used, 
as well as when the initial condition consists of two 
completely coherent parts. These dynamical properties 
indicate that the flow within a certain volume of the 
phase space is directed towards the amplitude chimera 
state. From the perspective of the amplitude chimera, 
there must exist some associated “stable directions”. 
However, even in the absence of any external perturba¬ 
tion, for all system sizes, the amplitude chimera states 
disappear after some time, and the system approaches a 
coherent oscillatory state. Accordingly, there must also 
exists at least one “unstable direction” in phase space. 
These findings can be explained by the structure of the 
phase space, which is schematically depicted in Fig. [15) 


The lifetime of amplitude chimeras in the deterministic 
case strongly depends on the choice of initial conditions 
as discussed in Sect. Ing In general the sensitivity of 
chimera states to the initial configurations is explained 
by the fact that classical chimera states typically coex¬ 
ist with the completely synchronized state, for which the 
basin of attraction is significantly larger. For amplitude 
chimeras, all our numerical results support the idea that 
amplitude chimera patterns can be seen as a saddle state 
composed of stable (solid lines in Fig. 15) and unstable 
(dashed lines Fig.[^ manifolds. The set of initial condi¬ 
tions leading to amplitude chimeras can be represented as 
a volume restricted in phase space (green disk in Fig. 15). 
The observed amplitude chimera corresponds to trajec¬ 
tories starting from this set and passing the saddle-point 
from the stable direction towards the unstable manifold. 
The lifetime of an amplitude chimera, therefore, depends 
on the chosen trajectory: the closer to the saddle-point it 
gets, the longer is the lifetime. In other words, the tran¬ 
sient time is determined by the time the system spends 
in the vicinity of the saddle-point where coherent and 
incoherent oscillating domains coexist before it escapes 
to the in-phase synchronized regime along the direction 
of the unstable manifold. Such a phase space scenario 
explains the sensitivity of transient times to initial con¬ 
ditions since they determine the particular path the sys¬ 
tem takes. With increasing network size the dimen¬ 
sion of the phase space of the network increases, and it 
becomes more likely that the distance of the initial con¬ 
ditions to the stable manifold is larger, which leads to 
decreasing transient times for passing the saddle-point. 
This explains the observed decrease of the lifetimes of 
the amplitude chimeras with N (FigjTT]). 




1-CD 

3-CD 

n-CD 

SYNC 


FIG. 16. (Color online) Map of dynamic regimes in the plane 
of coupling strength a and coupling range P for noise inten¬ 
sities: (a) D = 0, (b) T) = 5 • 10“^^, (c) D = b ■ 10“^^, (d) 
D = 5 • 10“^. Color code: 1-cluster chimera death (1-CD), 
3-cluster chimera death (3-CD), multi-cluster chimera death 
(n-CD, n>3), in-phase synchronized oscillations and coher¬ 
ent traveling waves (SYNC). Initial condition: snapshot of an 
amplitude chimera calculated for T) = 0, P = 4, cr= 14, t=150. 
Maximum simulation time: t = 5000. Parameters: N = 100, 
A = l, uj = 2. 


Our numerical investigations of the stochastic model 
Eq. <§ show that Gaussian white noise dramatically re¬ 
duces the impact of initial conditions on the lifetime of 
amplitude chimeras. In more detail, we have tested a set 
of realizations of initial conditions which lead to signifi¬ 
cantly different lifetimes of amplitude chimeras without 
random forcing. In the presence of relatively weak noise 
D = b- 10“^^ all realizations result in amplitude chimeras 
with similar lifetime. This again supports our view of the 
amplitude chimera as a saddle-point, and allows for the 
following explanation. The stochastic force, which con¬ 
tinuously perturbs the system, makes it randomly switch 
between different trajectories close to the saddle-point. 
Therefore, the system’s dynamics is not determined by a 
single trajectory anymore, but rather affected by a set 
of trajectories belonging to the Wdimensional hyper¬ 
sphere. This reduces the sensitivity of the amplitude 
chimera lifetime to specific initial conditions. In Fig.p!5| 
the impact of noise is illustrated by yellow shading, de¬ 
noting the stochastic forces applied to the system at one 
instant of time. 


V. MAPS OF DYNAMIC REGIMES 

For a large range of the coupling parameters a and 
P we calculate the asymptotically stable state and the 
transient time of amplitude chimeras for A = 100. For 
each choice of (a, P) we start with the same amplitude 
chimera configuration as initial condition. For an 
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FIG. 17. (Color online) Boundary between the oscillatory 
regime and the chimera death regime for different noise inten¬ 
sities D, extracted from the maps of dynamic regimes shown 
in Fig.[^ 


exemplary initial condition, the results belonging to four 
different noise intensities are shown in Fig.[^ For very 
small coupling strength a or very small coupling range 
P the asymptotic states are coherent oscillatory states, 
either in-phase synchronized oscillations or traveling 
waves (dark blue region, labeled SYNC). For very 
small coupling range, we observe amplitude chimeras 
as transients. For larger a and P, we find chimera 
death states (yellow, orange, and blue regions) with 
one coherent domain (1-CD), or for slightly smaller P, 
with three (3-CD), or more (n-CD with n > 3) coherent 
domains. For all noise intensities, there exists a chimera 
death regime (1-CD, 3-CD, n-CD), as well as a coherent 
oscillatory regime (SYNC). 

The regime of chimera death states is characterized 
by high multistability, characterized by different initial 
conditions. The boundary between the oscillatory regime 
and the chimera death regime is roughly independent 
of the particular amplitude chimera snapshot used as 
initial condition. In contrast, for many values of (cr, P), 
the particular type of chimera death depends on the 
realization of the initial condition. Note that there is 
nevertheless a clear tendency that the m-CD patterns 
with {m < k) are generally found for larger coupling 
ranges than the k-CD states {k,m G {l,3,n}). This 
tendency is especially pronounced for large coupling 
strengths. Noise influences the dynamic regimes in 
different ways. First, the boundaries between the 
different cluster types of chimera death appear to be 
almost unaffected by the applied external noise. We 
do not observe any noise-induced switching between 
the different types of chimera death. The applied 
noise does not infiuence the asymptotic chimera death 
state. Second, with increasing noise intensity, the 
boundary between the oscillatory regime and the oscil¬ 
lation death regime is shifted towards higher coupling 
strengths. This means that the stochastic force pushes 



FIG. 18. (Golor online) Transient times of amplitude chimeras 
ttr in the plane of coupling strength a and coupling range 
P, for the noise intensities: (a) P = 0, (b) P = 5 • 10“^^, 
(c) P = 5 • 10“^^. System parameters, initial condition and 
simulation time as in Fig. |16| 


the system out of the deterministic inhomogeneous 
steady state into the basin of attraction of the stable 
coherent oscillatory state, and induces oscillations in 
a parameter regime where in the absence of noise the 
steady state is a stable asymptotic solution. The size 
of this parameter regime depends on the applied noise 
intensity. In order to facilitate the comparison, the 
boundaries between the oscillatory regime and the 
chimera death regime are depicted for different noise 
intensities in Fig.[^ Generally, the stronger the applied 
noise is, the smaller is the regime of chimera death states. 


In the oscillatory regime, we observe transient ampli¬ 
tude chimeras. In Fig.[^ their lifetimes are depicted, 
obtained from the same simulations described above. 
One can see that generally the transient time decreases 
with decreasing coupling strength, and with increasing 
noise intensity, as shown already in Figs.[^ and M for 
a restricted range of coupling parameters. Note that 
Fig.f^b) (P = 5 • 10-11) and (c) (P = 5 • IQ-^i) look 


very similar up to rescaling of the transient times. This 
illustrates that the impact of the applied noise upon 
the dynamics is rather independent of the strength and 
range of the coupling. 


In the deterministic case (Fig.j^a)) there is a regime 
of high values of the coupling strength, at the border 
between the oscillatory regime and the chimera death 
regime, where the transient amplitude chimeras last 
longer than the maximum simulation time of t = 5000 
(bright orange). For several values of (cr,P) in this 
region, we have simulated much longer time series until 
t = 40000 (more than 12700 oscillation periods T), 
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and have found that the amplitude chimeras persist. 
However, they disappear much earlier as soon as a tiny 
amount of external noise is applied. This indicates 
that the amplitude chimera states are still unstable in 
this region. The extremely long transient times might 
simply be related to our choice of initial conditions in 
the deterministic system. 


VI. CONCLUSION 

In this work, the robustness of chimera states with 
respect to noise has been investigated for a paradigmatic 
network of oscillators. We have presented numerical 
results demonstrating that transient amplitude chimeras 
and chimera death states in a ring network of identical 
Stuart-Landau oscillators with symmetry-breaking 
coupling continue to exist in the presence of Gaussian 
white noise. 

Transient amplitude chimeras occur in the same range 
of coupling parameters as in the deterministic case. The 
key quantity we use to measure their robustness is the 
transient time. The latter decreases logarithmically 
with the applied noise intensity. For a constant noise 
intensity, the transient times increase with the coupling 
strength up to a saturation value. The width of the inco¬ 
herent domains relative to the overall system size is not 
affected by external noise, and is independent of the used 
initial condition, and the total system size. It increases 
linearly with the coupling strength as well as with the 
coupling range. The amplitude chimeras appear as long 
lasting transients in systems of all sizes, in contrast to 
classical phase chimeras whose transient times increase 
exponentially with the network size IHIIT]. The ampli¬ 
tude chimera lifetime decreases with the network size, 
but approaches a finite value in the thermodynamic limit. 

The amplitude chimera lifetimes depend sensitively 
on the particular realization of the randomized initial 
condition. We have introduced a class of specially 


prepared random initial conditions that produce long 
lasting amplitude chimeras. We have shown that initial 
configurations that fulfill a symmetry conditions which 
is also found in oscillation death patterns result in the 
longest living amplitude chimera transients. 

We have further introduced of set of global order 
parameters which are capable to detect the transitions 
from transient amplitude chimeras to completely in- 
phase synchronized oscillations or coherent traveling 
waves. Moreover, they can be used as a measure of the 
width of the incoherent domains of amplitude chimeras. 

The chimera death patterns also persist under the 
impact of stochastic forces. However, the coupling 
parameter regime where they occur is reduced with 
increasing noise intensity. The boundary between the 
coherent oscillatory regime and the chimera death 
regime is shifted towards higher values of the coupling 
strength. That means that the system favors oscillatory 
behavior for a larger coupling parameter regime. In 
contrast, this boundary appears to be independent on 
the particular realization of the initial condition. The 
number of clusters within the coherent domains appears 
to be unaffected by the external noise, but depends on 
the particular initial condition. 

Our numerical findings can be explained by the 
underlying phase space structure. More specifically, 
we propose that amplitude chimera states can be 
represented by saddle states in the phase space of the 
network. This elucidates the behavior of their lifetime, 
and explains that generally the initial conditions become 
less important under the influence of noise. 
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